Tradeoff between energy and error in the discrimination of quantum-optical devices 
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We address the problem ol energy-error tradeoff in the discrimination between two linear passive 
quantum optical devices with a single use. We provide an analytical derivation of the optimal 
strategy for beamsplitters and an iterative algorithm converging to the optimum in the general case. 
We then compare the optimal strategy with a simpler strategy using coherent input states and 
homodyne detection. It turns out that the former requires much less energy in order to achieve the 
same performances. 

PACS numbers: 



I. INTRODUCTION 



There are many contexts in which one wishes to use 
as little energy as possible in order to perform a desired 
task. Consider the case in which one wants to gather in- 
formation about an unknown device. The most general 
approach to the problem is to probe the target and to 
measure its response. This implies an undesired pertur- 
bation of the system. Often such perturbation depends 
on the energy of the probe, and one can ask what is 
the minimum amount of energy the probe must have in 
order to perform the desired task. Exploiting the distinc- 
tive features of quantum theory, such as entanglement, 
one could in principle achieve better performances with 
respect to a classical strategy with the same energy. For 
example, in the context of quantum illumination [H-Qj 
one discovers that the use of entangled light allows to 
enhance the detection and the imaging of an object. 

The task we are considering in this paper is the dis- 
crimination between two quantum optical devices (such 
as beamsplitters or phase shifters). An optical device is 
represented by a unitary transformation U acting on a 
system described by a Fock space H. The general prob- 
lem in which one wants to discriminate between two uni- 
taries has been widely investigated in the literature 
and can be summarized as follows. Suppose we are pro- 
vided with a single use Q of an unknown device randomly 
chosen in a set of two Ui and U2, with equal prior proba- 
bility. Our task is to distinguish between Ui and U2 with 
a probability of error not greater than a given threshold 
q. One sends an input state |'0) that undergoes the uni- 
tary evolution described either by Ui or U2- Finally, one 
applies the optimal POVM [l3| in order to distinguish 
the two possible output states Ui {ip) ot U2 \'ip) with min- 
imum error probability. This error probability depends 
on the input 1-0), so one can optimize on In the case 
at hand, Ui and U2 describe optical devices and the input 
state \tP) is a quantum state of light, with a well-defined 
mean value of the energy. Since in general many states 
allow to discriminate with an error probability Pe < q, 
one can wonder which state \^p*) accomplishes this task 
with minimum mean value of the energy. 



Here we analyze the energy-error tradeoff in the case 
in which the optical devices to be discriminated are lin- 
ear, passive and lossless. This discrimination is useful for 
example when reading classical digital information en- 
coded in the reflectivity of a media, such as conventional 
CDs or DVDs. This particular application has flrst been 
suggested in [ll| for the discrimination of quantum chan- 
nels affected by loss and noise. In particular, there it is 
shown that, for fixed mean number of photons irradiated, 
non-classical light can outperform any classical source in 
terms of the amount of information retrieved. Recently, 
the discrimination of lossy beamsplitters has been con- 
sidered in [12|. in the scenario where one has access to 
only one of the input and output modes and to a re- 
stricted class of input states. A similar scenario, namely 
the discrimination of a lossy quantum channel from an 
ideal one, has been studied also in [Tsj . 

The paper is organized as follows. In Section |ll] we 
review some basic notions of linear optics and introduce 
the problem. The analytical derivation of the energy- 
error tradeoff when the devices are beamsplitters is given 
in Section IIII Al while a numerical algorithm for the so- 
lution in the general case is presented in Section IIIIBI 
In Section IIII CI we analyze the energy-error tradeoff in 
the restricted scenario in which only coherent states and 
homodyne detections are available. We will then quan- 
tify the advantages that one has by adopting the optimal 
quantum strategy. Section IIVI concludes the paper with 
a discussion of the results. 



II. DISCRIMINATION OF LINEAR PASSIVE 
OPTICAL DEVICES 

A Af -modes quantum optical device [3, [l^ is de- 
scribed by a unitary operator U relating M input optical 
modes with annihilation operators on Hi, to M out- 
put optical modes with annihilation operators a[ on "Hii , 
where Hi denotes the Fock space of the optical mode i. 
We denote a.s V. — 0^ Hi. 

An optical device is called linear if the operators of 
the output modes are related to the operator of the input 
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modes by a linear transformation, namely 

N N 

a' := UaiU^ = ^ Ai^aj + ^ Sij-a], i = 1, . . . , iV. 

The above equation can be rewritten in the more compact 
form 



a' ' / \ a 
where S' is the 2N x 2iV scattering matrix defined as 

A B 



S 



B A 



(1) 



(2) 



[X is the complex conjugate of X), a = [ai^ . . .a^) is 
the vector of annihilation operators of the input mode, 
and analogously a' for the output modes. If S = in 
Eq. ^ the device is called passive and conserves the 
total number of photons, that is 



{'ip\N\^) = (v-ic/^ivt/iv-) 



(3) 



with N := a\ai the number operator on %. 

Suppose now that we want to discriminate between 
two linear optical passive devices Ui and C/2- If a single 
use of the unknown device is available, the most general 
strategy consists of i) preparing a bipartite input state 
p e B{T-L (E) JC) {JC is an ancillary Fock space with mode 
operators bi), ii) applying locally the unknown device and 
iii) performing a bipartite POVM 11 — {111,112} on the 
output state {Ux ®Ik:)p ~ {Ux ® Ik)p{UI ® Ik) [x can 
be either 1 or 2) 



Ux 



n 



(4) 



When the device is randomly chosen from the set 
{f/i,t/2} with equal probabilities pi — P2 — 1/2, the 
minimum probability of error in the discrimination can 
be proved to be [IB] 

Pe{p,Ui,U2) = ||(Wl®Z;c)p-(^2®I;c)p||l)• 
(5) 

where ||X||i = Tr[VxVC]. If we define as N/c = bjbi 
the number operator on the ancillary modes 6^, the en- 
ergy of the state p is proportional to 



E{p) 



-Tr[piN(g,lK.+I(g>NK:)], 



(6) 



while clearly the energy that flows through the unknown 
device is Ed{p) ■= ^ + Tr[p{N ® I^)]. 

Since we have Pe{{Ul ®Ik)p, I, U2UI) = Pe{p, Ui, U2) 
and Eq. implies E{{ul (S)Ijc)p) = E{p), we can 

restrict our analysis to the case in which Ui — I and 



U2 = U. We consider now the problem to find the mini- 
mum energy input state p* that allows us to discriminate 
between / and U with probability of error not greater 
than a given threshold q, that is 



p* = arg min E(p). 
^ps.t.PM<q 



(7) 



Sometimes one is more interested in minimizing the en- 
ergy Ej,{p) flowing through the device, rather than the 
total energy E{p) of the input state. In the following we 
show that no ancillary modes are required for optimal 
discrimination, so these minimization problems lead to 
the same optimal input state. 

The following lemmas allow us to simplify Eq. ([T]). 

Lemma 1. Without loss of generality the minimization 
in Eq. ([7]) can be rewritten as 



P 



arg mm 

p s.t. P^(p,U)^ 



E{p). 



(8) 



Proof. Suppose that for the optimal state p* one has 
Pe{p) < q- Since Pe{p) is a continuous function there 
exists a < a < 1 such that Pe((l — ct)p + oi |0) (0|) = q, 
where |0) (0| denotes the vacuum state. By observing 
that E{{l-a)p+a |0) (0|) < E{p) we have the thesis. □ 

Lemma 2. The optimal state achieving the minimum in 
Eq. ([7]) can be chosen pure. 

Proof. First, let us prove that the minimization of E{p) 
for a given value of Pe is equivalent to the minimization 
of the convex combination 



Fip)^pP,ip) + {l~p)Eip) 



(9) 



for a fixed value of p. Suppose that we find /5*that min- 
imizes E{p). It is then clear that for q := Pe{p*), E{p) 
gives the minimum possible value for the energy because 
any lower energy would decrease F . From Eq. ([5]) and 
the convexity of the trace distance, it follows that Pe{p) 
is a concave function of p. Since Pe{p) is concave and 
E{p) is linear, F{p) is a concave function of p and its 
minimum is attained on the boundary of its dominion, i. 
e. for pure states. □ 



Lemma [2] allows us to rewrite Eq. ^ as 

Pe. = ^ (1 - v^^\mu®hm?) 



(10) 



We now prove that no ancillary modes are required. 

Lemma 3. Without loss of generality, the minimum in 
([7]) is achieved without using ancillary modes. 

Proof. Due to Lemma [2] the input state can be written 
as IV') — I*) \Xi) where \i) is an orthonormal basis 

in H and \xi) are normalized states in JC. If we define 
lip') :— Ci \i) |0) (|0) is the vacuum state), it is easy to 
verify that Pe(V') = Pe{^') while E{i:) > E{'ip'). □ 
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After these consideration and disregarding some irrel- 
evant constant factors, Eq. ([7]) can be rewritten as 



arg min {i>\N\'di), (11) 

l^>s.t.|(vi;7|^)l=K' 



where from Eq. (dni) it follows that K = ^jAq{l - q). 



III. THE OPTIMAL TRADEOFF 

A. The case of beamsplitters 

In this section we derive an analytical expression for 
the optimal energy-error tradeoff for the case in which 
the devices to be discriminated are two beamsplitters. A 
beamsplitter is a two-mode linear passive quantum opti- 
cal device whose scattering matrix S has the form 



A 
A 



A^SU{2). 



(12) 



In the following we will use the basis {|nm)} with respect 
to which U is diagonal i.e. 



A = 



e- 







-iS 



< 5 < vr 



(13) 



With this choice, for any lip) = J2'^,m=o'^nm \n,m), we 

nm) and it is easy to observe 
a„™|V*("-")and(V'|A^|V') = 



have U \nm) 



that(V|f/|V')=E^ 

E 



n,m=0 l^nml 



'{n + m). We notice that both these ex- 
pressions only depends on the squared modulus of the 
coefficients a„m and so we can assume anm to be real 
and positive. The assumption that the devices are beam- 
splitters allows us to simplify the structure of the opti- 
mal input state for the tradeoff. First, we show that it is 
not restrictive to consider superpositions of the so called 
NOON [i3 states. 



Lemma 4. If U is a beamsplitter, the optimal state \^* 
in Eq. (jlip can be taken of the form 



oo 

E< 

n=0 



I0n) , !</'«) 



i(|n,0) + |0,n)) (14) 



Proof. We show that for any state 



n.m=0 



\n,m) there exists a state \^') 



EZo^^il^i) such that iij'lNli;') < {^\N\^) 
and \{i;'\U\i(;')\ < \{i;\U\i^)\. Upon defining 



= E|„-ml=i l"nm|^, One Can verify 



{4''\Nm= alJn-m\<{^\N\ij), (15) 

n,m— 



"mnCOs((S|n-TO|) 



< 



<\WUW 

that proves the statement. 



(16) 



□ 



From Eq. p4)) it follows that the expectation value 
of U over \ip) is real. Following an argument similar 
to the one we used to prove Lemma [T] , the constraint 
I (V'l 1^) I = K can be changed into {ip] U\tp) — K and 
Eq. (fTTj) becomes 

W) = ais min (^|iV|-0). (17) 

\^)s.t.(-4,\U\^li)=K 

One can observe that for any state j-i/') of the form of Eq. 
(fT4|) . the constraint "-"^n cos((5n) = K implies that 

there must exist at least one non-null coefficient Q!„ for n 
s.t. cos{5n) < K and at least one for n s.t. cos((5n) > K . 
In the following lemma we will prove that not more than 
two non-null coefficient q;„ are needed. 

Lemma 5. The optimal state that minimize Eq. 

(|17p can be taken of the form 



W) 



(18) 



Proof. Let us consider the optimal state j?/;) — 
with (^|C/|^) = K and {iP\N\4,) = N^i^. 
Suppose now that the set a :— {an} has TV > 3 elements. 
Then there must exist ni and n2 such that a„j , a„2 
and cos{Sni) < K < cos{Sn2). It is then possible to de- 
fine \x) ■■= Pn^ \(t>m) + /5„2 I0n2> such that (xl U \x) = K. 
Furthermore, we can define |^) := (1 — e)^^/^ J2n 7" 1^")' 
where 



7n 



a„ if n 7^ ni, n2 



and e < min(a„j//3„j^, a„2//3n2)- We notice that 
(fl U \0 = K, and 

N^in = e{x\N\x) + il-e){^\N\0. (19) 

If (x|A^|x) = Nn^in the statement follows with = 
Ix). If {x\N\x) ^ iVmin, either (xl^lx) < ^mi„ or 
(CI ^ 10 < -^min, that contradicts the hypothesis that 
lip) is the optimal state. □ 

Since we both require (-01 [/ = if and a^^ +«n2 = 1 
the expression of the coefficients , and of the mean 
value ("01 N \tp) are fixed by the choices of ni, n2 and K, 
i.e. 



cos{Sn2) - K _ . 2 



cos{6n2) — cos(fei) 



(20) 



(0|JV|0) 



n2 cos((5ni) — ni cos{Sn2) + K{ni — 712) 



cos((5ni) — cos((5n2) 



(21) 



It is now convenient to rephrase the problem at hand 
in a geometrical way (see Fig. [1]). Let us introduce the 
map f : N — IR2, defined as f(n) — {cos{Sn),n). We 
can associate to any couple ni and n2 the line £„i,„2 •= 
{tf{ni) -I- (1 — t)f{n2)\t £ R} and the segment lni,n2 ■= 
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{tf{ni) + (1 - t)f{n2)\0 < t < 1}. It is easy to prove 
that the expectation value N\ip) in Eq. ([^0]) can be 
rewritten as 



ni ,712 



(22) 



where we defined the line := {{K,y)\y e R} and the 
mapping (a, = 6. For any K it is then possible to 
define a partial ordering among the segments lni.n2 
for which cos{Sni) < K < cos{Sn2), as follows: 

^ni,n2 ^n[^n2 if (^ni,Ti2 1^ '^K^y ^ {^n{,n2 ^ ^K^y (23) 
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FIG. 1: (Color online) Geometrical representation of the op- 
timization problem. For fixed values ni and n2, each state 
of the form (|18p is represented by a point of the line segment 
lni,n.2 which extrema are (cos(5ni), ni) and {cos{5n2) , n2) ■ 
The X and y coordinates represent the expectation values 
{tp\U\^p) and {%l)\N\tp), respectively. From the picture it is 
clear that the line segment /o,n* lies below any other line seg- 
ment Zni,n2, and so it identifies the family of optimal states, 
as we prove in Prop. [1] 



nonempty. Consider now the lines Ln'.o and Ln'^m- It 
is easy to verify that for all m we have {Ln'^m ^fi)y > 
{Ln',o(^ri)y which implies {Ln',m(^rK)y > (-Z^n'.o H r/^)y 
for all cos((5n') < K < I and finally In'.m >k In'fi 
for all cos{Sn') < K < cos((5m). Similarly one can 
verify the bound ln',o <k 1-^/5,0 <k lni.n2 that holds 
for cos((5ni) < K < cos{Sn2), n' € O. We can then 
restrict the optimization over the finite set of states 
S := I lip) = an' + ao |00) ,n' e Q}. As a con- 
sequence, the mean value {tp\ N j-f/;) becomes 



1 — cos{Sn) 



(25) 



The right hand side of Eq. ((25|) can be proven to be 
a convex function [18] for tt/2 < 6n < tt and achieves 
its minimum for h minimum positive solution of 6n = 
tan((5ri/2). Since n is in general not integer the optimal 
value n* is given by argmin^^^j jf^^ N |V'*)- CH 

Figure [2] shows the optimal energy-error tradeoff ob- 
tained with the discrimination strategy of Prop. [TJ 



3.5 



5=71/6 (upper line) 
6=jt/4 (middle line) 
6=n/3 (lower line) 




0.5 



FIG. 2: (Color online) Optimal tradeoff between the energy 
E and the probability of error Pe in the discrimination of I 
and U = exp(i(5a|ai — Sala2)), for various values of S. 



We are now ready to proof the following Proposition. 
Here \_x\ ( \x~\ ) denotes the minimum (maximum) integer 
number greater (smaller) than x. 

Proposition 1. For any K, < K < 1, the optimal 
state lip*) that minimizes Eq. (|17p is 

|V'*) = a«*|0n.)+«o|OO), (24) 

where n* — argmin^fjj p^^^ {'>p*\ N [ip*) , with h the mini- 
mum positive solution of Sn = tan((5n/2). 

Proof. First let us introduce the set := {n € N|7r/25 < 
n' < tt/S}. Since S < tt we observe that ft must be 



B. An iterative algorithm for the general case 

Here we provide an iterative steepest-descent algo- 
rithm [l^ to find a state \ip) that solves the optimization 
problem in Eq. . With the same argument we use in 
the proof of Lemma [2] we can rephrase the optimization 
problem in Eq. (|lip as the minimization of the convex 
combination 

C(V') := p (VI N IV) + {l-p)\ (VI U IV) p. (26) 
We are now ready to introduce the iterative procedure. 
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Proposition 2. The following algorithm converges to a 
state \ip) that is optimal according to Eq. (|26p. Take an 
arbitrary state jV'")- Given IV""); evaluate by the 

following steps: 

1. Evaluate the derivative of the figure of merit C 



[pN+il-p)i{r\U\r)U^+ (27) 
+ (V'"|J7^|V'")C/)]|V'")- (28) 



2. Pick up a small enough positive a and evaluate 

5C(V'") 



n+l 



3. Normalize 



(1 - a) IV'") - a- 



(29) 



n+l 



according to 



|^"+i) = |^"+i|-i 



n+l 



(30) 



C. Discrimination of passive devices with limited 
resources 

Here we consider the minimum energy discrimination 
tliat makes use of colierent input states \ai) and homo- 
dyne detections X^p. to discriminate a single use of a n- 
modes passive linear optical device randomly chosen in 
the set {/, U} with equal probabilities. 



















Xip2 













(31) 



If we consider coherent input states \ai) on mode i the 
global input state is |^) = 0^ \ai) which corresponds to 
an energy value 

E{0 (ei E\0 = l + {^\N\0 = l+J2 (32) 



Proof. The algorithm in Prop. [2] is a steepest-descent 
algorithm: we move the state in the direction of the gra- 
dient of the figure of merit C{tp), so by construction one 
has C(V'"+^) < C(V'"). □ 

The parameter a controls the length of each iterative 
step, so for a too large an overshooting can occur. This 
can be kept under control by evaluating the figure of 
merit C{i>) at the end of each step: if C('0) increases 
instead of decreasing, we are warned that we have taken 
a too large. Figure [3] shows the optimal energy-error 
tradeoff obtained with the algorithm of Prop. [5] 



5i=7r/4, b2=-nlA (upper line) 
6^=7r/2, 82=0 (lower line) 




Since for any passive linear device V we have that 
= ®i ^ith coherent state, we can 

assume U to be diagonal, i.e. U = e'^'^^"'. The evo- 
lution of 1^) under the action of U is then given by 



u\0 



(33) 



A quantum homodyne detection X^ is described jlj, 
[13 by the POVM {|a;, ip) (x, ip\} where |a;, ip) are the 



eigenstate of the quadrature e 



The prob- 



ability of outcome x when the system is prepared in a 
coherent state |a), a = e**^" \a\ is given by the Gaussian 



p^{x\a) = \{a\x,ip)[' 



1 



/27r 



(34) 



We notice that p^{x\a) depends on Lp only trough the 
sum Lp + (pa- We can then fix = and vary only the 
Q,;. conditional probabilities of outcome x = {xi) given / 
or U are n-dimensional Gaussians, namely 

p(x|/) = (27r)-"/2elx-volV2^ 



0.5 



p(x|[/) = (27r)-"/2g|x-vi| /2^ 

vo = (2Reai), vi = (2 Re e'^'a,), 



(35) 



FIG. 3: (Color online) Optimal tradeoff between the energy 
E and the probability of error Pe in the discrimination of two 
linear passive optical devices I and U . We consider the case 
in which i7 is a two mode device, namely U = exp(i((5ia|ai + 
52a\a2)). 



Any classical postprocessing of the outcome x can be 
described by a function g(X|x) that evaluates to 1 if one 
guess the unitary X from outcome x, and to otherwise, 
with X — I,U. The probability of error is given by 



dxp(x|/)g(;7|x) +p(x|C/)g(/|x), (36) 
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and thus the optimal postprocessing is 

n(XM-l 1 ifMx|X)>p(x|r) 
<1[^W - i ifp(x|X) < p(x|X) 



6=ji/4 (a) 



(37) 



Inserting Eq. (|37l) and Eq. ([351) into the expression p6| , 
the probabihty of error becomes 



(27r) 



-n/2 



where we defined the set 

A = {xs.t. |x- vop > |x- vip}. 



(38) 



(39) 



Within this framework it is more convenient to fix the 
amount of energy, that is the average number of pho- 
tons 77, and find the input state |^*) that minimizes the 
probabihty of error in the discrimination, i.e. 



ir) = arg min P^{^). 



(40) 



With a httle machinery it is possible to prove that 
Pe(l^)) is a non-increasing function of |vo — vip and 
then the minimization of Pe(|^)) can be rephrased as a 
maximization of |vo — Vip. We have then 

|vo-vi|2 ^4^[Re(aO-Re(e'^aO]' = 

i 

= iJ2^{cos{<j>i) - cos{(j)i + S,))\a^\f < 

i 

< 4^[2sin((5i/2)|ai|]2 < 16 sin^ ((5/2)77 

(41) 

where 6* := argmax^. \Si\, and i* labels the correspond- 
ing mode. The bounds in Eq. ([1T|) are achieved for 



ID = (g) |0,) ® \a* 



(42) 



y^exp(7 ^ ). The corresponding optimal 



where a* 
discrimination strategy is 



a*. 








Xo 


cm— 













(43) 



where — fT) means that the corresponding mode is dis- 
carded. With this choice of the input state the probabil- 
ity of error becomes 



(44) 



where $(a;) — dtexp^—t"^ /2) denotes the normal 

cumulative distribution function. 



Coherent states - tiomodyne detections strategy (upper line) 
Optimal strategy (lower line) 




FIG. 4: (Color online) Optimal tradeoff between the energy 
E and the probabihty of error Pe in the discrimination of I 
and U = exp{i{Sa\ai — Sa2a2)) {S = 7r/4 in (a) and S = 7r/12 
in (b)). The upper line represents the discrimination with 
coherent states and homodyne detections, while the lower line 
represents the optimal discrimination. Comparing (a) and 
(b), we notice that the improvement provided by the optimal 
strategy increases as S decreases. 



From Eq. (|44)) one can obtain the tradeoff between the 
energy and the probability of error, which is plotted in 
Fig. m for some choices of Ui and 1/2- If we consider the 
case in which we want to discriminate a 50/50 beamsplit- 
ter from the identity, one can notice that, for Pg = 0.1, 
the coherent state - homodyne detection discrimination 
strategy requires a factor of ^ 4 more photons that the 
optimal strategy. Moreover, this factor increases as the 
two devices get closer, i. e. for small values of 6. For ex- 
ample, when 5 = 7r/12, the factor is ~ 12. As expected, 
one notice that this factor increases when the probability 
of error decreases. 
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IV. CONCLUSION 

In this paper we studied the energy-error tradeoff in 
the discrimination of two linear passive optical devices. 
We derived the optimal strategy when no restrictions 
in the input state and in the final measurement are as- 
sumed. It is shown that the input state can be taken as 
pure and no ancillary modes are needed. For the beam- 
splitter case the input state is proved to be a coherent 
superposition of a NOON state and the vacuum 

1 00). It is worth noting that the choice of | depends 
only upon the reflectivity of the beamsplitter. We pro- 
vided an iterative algorithm to solve the problem in the 
general case. 

In Section llll Cl we considered the practical scenario in 
which one performs the discrimination with coherent in- 



put states and homodyne detections. This strategy, when 
compared with the optimal one, turns out to be largely 
suboptimal. For example, for a 50/50 beamsplitter the 
energy required for the discrimination when the proba- 
bility of error is 0.1 is 4 times smaller with the optimal 
strategy than with the coherent-homodyne one. 
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